Quantum logical gates with four-level SQUIDs coupled to a superconducting resonator 
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We propose a way for realizing a two-qubit controlled phase gate with superconducting quantum 
interference devices (SQUIDs) coupled to a superconducting resonator. In this proposal, the two 
lowest levels of each SQUID serve as the logical states and two intermediate levels of each SQUID 
are used for the gate realization. We show that neither adjustment of SQUID level spacings during 
^ . the gate operation nor uniformity in SQUID parameters is required by this proposal. In addition, 

O ' this proposal does not require the adiabatic passage or a second-order detuning and thus the gate 

' is much faster. 

00 '. PACS numbers: 03.67.Lx, 42.50.Dv, 85.25.Dq 

' Introduction. — Superconducting devices including cooper pair boxes, Josephson junctions, and superconducting 
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quantum interference devices (SQUIDs) have appeared to be among the most promising candidates for scalable 
quantum computing, due to design flexibility, large scale integration, and compatibility to conventional electronics 
[1-3]. In the past few years, for SQUID systems, many theoretical methods for realizing a single-qubit gate and a 
' two-qubit controlled-phase (CP) or controlled-NOT (CNOT) gate have been presented [4-15]. 

For realizing a two-qubit CP gate with SQUIDs, several methods have been proposed based on cavity QED technique 
[8-15]. These methods are of importance in building quantum logic gates and open a new avenue for the physical 
realization of quantum information processing with SQUIDs in cavity QED. However, we note that these methods 
^ have some disadvantages. For instances: (i) the methods presented in [8,9] require adjustment of the level spacings of 
^vq SQUIDs during the gate operation, thus decoherence caused due to the adjustment of level spacings may pose a severe 
i problem; (ii) the methods proposed in [10,11] require slowly changing the Rabi frequencies to satisfy the adiabatic 
CO ' passage and the approaches introduced in [12-14] require a second-order detuning to achieve an ofF-resonant Raman 
[ coupling between two relevant levels; note that when the adiabatic passage or a second-order detuning is applied, 
psj ■ the gate becomes slow (the gate time is on the order of one microsecond to a few microseconds [11,13]); and (iii) the 
' proposal reported in [15] employs a two-mode resonator/cavity as well as a second-order detuning between the two 
> cavity modes; technically speaking, the requirement for a SQUID interacting with more than one cavity or resonator 
' modes is difficult to meet. In addition, it is noted that although two-qubit CNOT, CP, or iSWAP gates have been 
J> , experimentally demonstrated in superconducting charge qubits, flux qubits, and phase qubits [16-18], to the best 
of our knowledge, no experimental demonstration of a two-qubit gate with SQUID qubits in cavity QED has been 
reported. 

In this paper, we present an alternative method for implementing a two-qubit CP gate with two SQUIDs coupled 
" ■ ' to a superconducting resonator. As shown below, this proposal has the following advantages: (a) there is no need for 
adjusting the level spacings of SQUIDs during the gate operation, thus decoherence caused by tuning the SQUID level 
spacings is avoided; (b) neither slowly changing the Rabi frequency nor the use of second-order detuning is required, 
thus the gate is significantly faster (as shown below, the operation time of the gate is on the order often nanoseconds); 
and (d) only one mode of the resonator is employed. In addition, this proposal does not require identical coupling 
constants of each SQUID with the resonator and thus is tolerable to inevitable nonuniformity in device parameters. 
We believe that this work is of interest because it avoids most of the problems existing in the previous proposals. 

Basic theory. — The SQUIDs considered throughout this paper are rf SQUIDs each consisting of a Josephson tunnel 
junction enclosed by a superconducting loop. The Hamiltonian for an rf SQUID, with junction capacitance C and 
loop inductance L, can be written in the usual form [19] 

(^-^.)^ /^^N 
2C 2L \ 

where the magnetic flux threading the ring, and Q, the total charge on the capacitor, are the conjugate variables 
of the system, is the static (or quasistatic) external magnetic flux applied to the ring, and Ej = Ic^o/2tt is the 
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FIG. 1: (a) SQUID-resonator resonant interaction, (b) SQUID-resonator off-resonant interaction. The difference between the 
level spacings in (a) and the level spacings in (b) can be achieved by choosing different device parameters for SQUIDs. 

Josephson coupling energy, where Ic is the critical current of the junction and $o = h/2e is the flux quantum. 

A) . SQUID-resonator resonant interaction. Consider a SQUID (say SQUID a) coupled to a single-mode resonator 
and driven by a classical microwave pulse. The SQUID is biased properly to have four lowest levels, which are denoted 
by |0) , |1) , |2) , and |3) , respectively [Fig. 1(a)]. The resonator mode is resonant with the |2) o |3) transition but 
decoupled (highly detuned) from the transition between any other two levels, which can be readily achieved by 
adjusting level spacings of the SQUID [8,20]. In the interaction picture, the interaction Hamiltonian for the SQUID 
and the resonator mode, after making the rotating-wave approximation, can be written as [8] 

= ft(5aC+|2)^(3|+H.c.). (2) 

Here, the subscript a represents SQUID a; c+ and c are the photon creation and annihilation operators of the resonator 
mode with frequency ojc] ga is the coupling constant between the resonator mode and the |2) o |3) transition of SQUID 
a. The initial state |3)^ |0)^ and [2)^^ of the system, under the Hamiltonian (2), evolve as follows 

|3)J0)^ ^ cosg,t|3)jO)^-zsin5,i|2)Jl)^, 

|2)J1)^ ^ cosg,t|2)Jl)^-zsin5,i|3)jO)^, (3) 

where |0)^ and are the vacuum state and the single-photon state of the resonator mode, respectively. 

B) . SQUID-resonator off resonant interaction. Consider a system composed of a SQUID (say SQUID 6) and a 
single-mode resonator. Suppose that the resonator mode is off-resonant with the |2) o |3) transition (i.e., Ac = 
^^32 — Wc ^ gb) while decoupled from the transition between any other two levels of SQUID h [Fig. 1(b)]. Here, Ac 
is the detuning between the |2) o |3) transition frequency a;32 of SQUID b and the resonator mode frequency Wc, 
and gb is the coupling constant between the resonator mode and the |2) |3) transition. The effective interaction 
Hamiltonian in the interaction picture can be written as [21,22] 

^e = fi£(|3),(3|-|2),(2|)c+c, (4) 

where the subscript b represents SQUID h. 

From the Hamiltonian (4) , it is straightforward to see that if the resonator mode is initially in a single-photon state 
1)^, the time evolution of the states of the system is then given by 

|2)Jl)c ^ e^^^/^M2)Jl)c, 

|3)Jl)c ^ e-^^*/''M3)Jl)c, (5) 

which introduces a phase shift e^Sb'^l^" to the state |2) while e"*^''*/^'^ to the state |3) of the SQUID, when the resonator 
mode is in the state |1)^ . Note that the states |2}^ |0}^ and |3)^ |0)^ remain unchanged under the Hamiltonian (4). 

In the following gate operations, we will need this resonant interaction between the pulse and SQUIDs. Note that 
the resonant interaction between the pulse and the SQUIDs can be completed within a very short time, by increasing 
the pulse Rabi frequency (i.e., by increasing the intensity /amplitude of the pulse). 

Two-qubit CP gate. — Let us consider two SQUIDs a and b. By choosing different device parameters for each SQUID, 
SQUIDs a and b can have the four- level configurations as depicted in Fig. 1(a) and Fig. 1(b), respectively. The two 
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FIG. 2: Illustration of SQUIDs interacting with the resonator mode and/or the microwave pulses during the gate performance. 
The figures on the left (right) side correspond to SQUID a (6). 



logic states of a SQUID qubit are represented by the two lowest levels |0) and |1) , while the two intermediate levels |2) 
and 1 3) of each SQUID arc utilized for the gate realization. For the notation convenience, wc here denote the ground 
state (the first excited state) as level |1) (|0)) for SQUID b [Fig. 2(a',b',c',d',e')]. We suppose that the resonator 
mode is resonant with |2) -f-)- |3) transition of SQUID a while off- resonant with the |2) ^ |3) transition of SQUID 
b, which can be reached by prior adjustment of the level spacings of SQUIDs a and b. In addition, wc assume that 
the resonator mode is initially in the vacuum state |0)^. The notations Wg^-*, W20) and UJ21 involved in the following 
gate operations are the |1) •(->■ |3) transition frequency, the |0) O |2) transition frequency, and the |1) O |2) transition 
frequency of SQUID i (i = a,b). 

The operations for realizing a two-qubit CP gate are listed as follows: 

Step (i): Apply a microwave pulse (with a frequency LJfj,^ = and a phase (j) = tt) to SQUID a for a time interval 

ti = tt/ {2^13) [Fig. 2(a)], to transform the state to i \3)^ . Then, wait for a time interval t[ = tt/ {2ga) during 
which the |2) o |3) transition of SQUID a resonantly interacts with the resonator mode [Fig. 2(b)], to transform the 
state 13)^ JO)^ to -i \2)^ , as shown in Eq. (3). 

It can be found that after this step, the following transformation is obtained: 



|l)jO)/nN|3)jO), ^"4*M2)J1),. (6) 

On the other hand, the state ]0)^ |0)^ remains unchanged. 

Step (ii): Apply a microwave pulse (with a frequency = uj^ and a phase (j) = to SQUID a [Fig. 2(c)] 
while cL microw&vB puls6 (with. & frcquGncy uj^w — ^21^ 

and a phase (j) = -tt/2) to SQUID b [Fig. 2(b')]. The Rabi 
frequency for the pulse applied to SQUID a is ^02 while the Rabi frequency of the pulse applied to SQUID b is fli2- 
We set r2o2 = ^^i2j which can be achieved by adjusting the intensities of the two pulses. After the pulse duration 
t2 = 7r/(2rio2) = 7r/(2f]i2), the state |2) (|0)) of SQUID a is transformed to the state JO) (- 12)) while the state |1) of 
SQUID b is transformed to the state ]2). 

Step (iii): Wait for a time ts. Note that in the case when the resonator mode is in the photon state ]1)^ , the levels 
|2) and 13) of SQUID a are not populated after the above operations. Therefore, there is no coupling between the 
resonator mode and SQUID a. The resonator mode is off-resonant with the 12) o |3) transition of SQUID b [Fig. 
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2(a')]. It can be seen from Eq. (5) that for ~ nAc/g"^, the state |2)(, changes to — |2)(, On the other hand, 
the state |0)j |0)^ , |0)(, , and |2)jj |0)^ remain unchanged. 

Step (iv): Apply a microwave pulse (with a frequency w^^, = uj!^ and a phase (j) = —Tr/2) to SQUID a [Fig. 2(c)] 

while a microwave pulse (with a frequency w^^j = 1^21 ^^'^ ^ phase cp = it/2) to SQUID b [Fig. 2(b')]. Like step (ii), 
we set r2o2 — ^12- After the pulse duration t2 given in step (ii), the state |0) (|2)) of SQUID a is transformed to the 
state |2) (— |0)) while the state |2) of SQUID b is transformed back to the state 

Step (v): Perform an inverse operation of step (i) [Fig. 2(a,b)]. That is, wait for a time interval t'l given in step 
(i), during which the |2) -o- |3) transition of SQUID a resonantly interacts with the resonator mode; and then apply 
a microwave pulse (with a frequency = oj^ and a phase </> = tt) to SQUID a for a time interval ti given in step 
(i). It can be verified that after this step, the following transformation is achieved: 

|2)Jl)/"^''^-.|3)jO)/n*Ml)jO),. (7) 

On the other hand, the state |0)^ |0)^ remains unchanged. 

The states of the whole system after each step of the above operations are summarized in the following table: 

|00)|0), |00)|0), |20)|0), |20)|0), 

|01)|0)^ Step(i) |01)|0)^ StepQi) |22)|0)^ StepQii) |22) |0)^ 

|10)|0), |20)|1), 100) 100) 

|11)|0), |21)|1), 102) |1)^ -|02)|1), 

|00) |0), |00) |0), 

StepQv) |01)|0)^ Stepfv) |01) |0)^ 

|20)[1), |10)|0), ' 

-|21)|1), -|11)|0), 



(8) 



where \kl) is abbreviation of the state \k)^ of SQUIDs (a, b) with k, I G {0, 1, 2}. It can be concluded from Eq. (8) 
that a two-qubit CP gate was achieved with two SQUIDs (i.e., the control SQUID a and the target SQUID b) after 
the above process. 

From the description above, it can be found that: (i) In contrast to the previous proposals [8,9], the method 
presented above docs not require adjustment of the level spacings of the SQUIDs during the gate operation; (ii) This 
method does not require slow variation of the Rabi frequency in contrast to [10,11]; (iv) Compared with the previous 
approaches [12-14], this method does not require a finite second-order detuning 6 = Ac — A^y, and thus the gate speed 
is improved by one order (here A^^u is the detuning of the pulse frequency with the transition frequency between 
the two associated levels of SQUIDs; for the details, see [12-14]); and (v) this method employs only one mode of the 
resonator, which is different from the previous proposal [15]. 

Discussion. — Let us give a brief estimate on the gate time. As shown above, the total operation time is 

T = 2ti + 2t[ + 2t2 +t3 = TT/ga + TrAc/fif^ + 7r/f2i3 + 7r/Oo2, (9) 

where O02 is equal to f2i2 (see steps (ii) and (iv) above). Without loss of generality, let us consider ^ 56 ^ 3.0 x 10^ 
s~^, which is available at present [9]. By choosing Ac = 10^6, O13 ~ Qq2 ~ lOflo) we have r ~ 12 ns. 
Several issues related to the gate operations above need to be addressed as follows: 

i) The level |3) of SQUID a is occupied in steps (i) and (v). Since only SQUID-pulse resonant interaction and 
SQUID-resonator resonant interaction are used in steps (i) and (v), the operation time ti +t[ in step (i) or (v), equal 
to tt/ {2Q13) + tt/ {2ga) , can be significantly shorten by increasing the pulse Rabi frequency O13 and the coupling 
constant ga- Alternatively, one can design the SQUID a to have a sufficiently long energy relaxation time 7^"'^ for the 
level |3) . By doing these, we can have 7^^ ^ ti + 1[, such that decoherence caused by the energy relaxation of the 
level 1 3) of SQUID a is negligibly small. 

ii) The occupation probability ps of the level |3) for SQUID b during step (iii) is given by [12] 

which need to be negligibly small in order to reduce the gate error. For the choice of Ac = Wgb, we have ps ^ 0.04, 
which can be further reduced by increasing the ratio of A^/gb- 

iii) For steps (i), (ii), (iv) and (v), the resonant interaction between the resonator mode and the |2) -f-^ |3) transition 
of SQUID a, involved during the application of the pulse, is unwanted. To minimize the effect of this unwanted 
interaction on the gate, the Rabi frequencies flis and f2o2 require to be much larger than the coupling constant ga, 
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FIG. 3: Average fidelity _F as a function of the Rabi frequency r2i2 (in unit of (?(,) for = lOpb. 



i.e., ^13, f2o2 ^ 9a- Note that this condition can be achieved by increasing J7i3 and r2o2 (i-e., via increasing the pulse 
intensity). 

iv) For either step (ii) or step (iv), when the SQUID h is in the state |2) and the resonator mode is in the single- 
photon state [1)^ , the unwanted off-resonant interaction between the resonator mode and the |2) o |3) transition of 
SQUID b induces a phase shift e^A^^I^" to the state |2) of SQUID 6, which will affect the desired gate performance. 
The effect of this unwanted SQUID-resonator off-resonant interaction on the gate can be made negligibly small as 
long as the condition 57i2 3> g^/Ac is met. In the following, we will give a discussion on the effect of this unwanted 
interaction on the fidelity of the gate. 

Suppose that the two SQUID qubits are initially in a generic state described by It/i (0)) — a |00) -f ji |01) + 7 |10) + 
9 |11}, where the coefficients satisfy the normalization. In the ideal case, it can be seen from Eq. (10) that after the 
five-step operations described above, the state IV" (0)) becomes j-^i^ (''")) = ol |00) |01) -I-7 |10) — 6* |11) . On the other 
hand, when the effect of the off- resonant interaction between the resonator mode and SQUID h is included during steps 
(ii) and (iv), one can easily work out the expression for the final state (r)) after performing the same operations 
above. To simplify our presentation, we will not give a complete expression for [■!/' (r)) due to its complexity. 

The fidelity is given by 



F = |(^.,(t)| V(t))|" 
= 1 - 2a; (1 +p- 



2 2 2\ 

q — r ) 



{l-q'-r'y + 2p^l + 



2 2\ 

q -r 



where 



16*1 , p = cos ((9, q : 



2^/^17 + ^4 



sm Lp, r 



simp, 



(11) 



(12) 



with s = gl/Ac and (p = Tr^fifj + s'^/i/ (2r2i2) . 

Eq. (11) shows that the fidelity is a function of x G [0, 1]. Thus, the average fidelity over all possible two-qubit 
initial states is given by 



F = 



F {x) dx 

[1+p'^ + q^ + + p^ (-1 + 2q^ - 2r^) + q^ {\ + 2r^)] . 



(13) 



It can be verified that when the unwanted "SQUID 6" -resonator off-resonant interaction in steps (ii) and (iv) is 
not considered (i.e., the case for = or s = 0), we have p = g = and r — 1, leading to F = 1 and F = 1. 

We have plotted the average fidelity F for the case Ac — lOgt (Fig. 3). One can see from Fig. 3 that the average 
fidelity F increases as the Rabi frequency i7i2 of the pulse applied to SQUID b becomes larger, and the F is ^ 1 when 
= O.Ggb- 

Conclusion. — We have presented a way to realize a two-qubit controlled phase gate with two SQUIDs, by the use of 
a microwave superconducting resonator. As shown above, in this proposal, (a) SQUIDs, which often have considerable 
parameter nonuniformity, can be used; (6) the adjustment of the level spacings, which is undesirable in experiment. 
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is avoided; and (c) neither the adiabatic passage nor a second-order detuning is needed and thus the gate can be 
performed much faster. 
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